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X x 
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M e 

P 

q 
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r t 

V 
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x,y,z 



a 

$ 

9 ,<$ 
n 



Drag Coefficient 

Lift differential coefficient 

Static moment differential coefficient 

Resultant force of perturbing influences 

Di ameter 

Force 

Acceleration of gravity 

Moment of inertia around lateral axis of projectile 
Moment of inertia around longitudinal axis of projectile 
Mass 
Moment 

Resultant moment of perturbing influences 
Angular velocity around Xg axis 
Angular velocity around y 0 axis 
Angular velocity around z 0 axis 

Distance between center of gravity and nozzle throat 

Ve locity 

Wind velocity 

Ground axis coordinates 

Angle of attack 

Side slip angle 

Angles defined in Chapt. 3 (2) 

Angles defined in Chapt. _ 3 (3) 

Angle rotated around Xg axis 



v 



M 



v 

IT 

P 

Suffix 

a 

B 

V 

x,Y, z 



Factor angle of trim of canard 
Trim angle of canard 

Factor angle caused by attitude of projectile 
Density of air 

Wind axis coordinate 
Body axis coordinate 
Velocity axis coordinate 
Rectangular axis 
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1. INTRODUCTION 



Solid body equations to calculate the trajectory of a projectile 
which is flying with thrust of ramjet under hypersonic Mach number are 
formulated. The projectile is gun launched. This projectile has an axi- 
symmetric configuration and is fin stabilized with four pop-out wings. 

The projectile has canards to control attitude. Ideally the projectile 
has no spin; however, due to effects from miscellaneous forces and moments 
which work on the body, it may acquire spin motion. 

In order to calculate the trajectory considering such conditions, 
solid body equations with 6 degrees of freedom are induced. One of the 
main problems which will occur during flight motion of the projectile is 
trajectory errors which are caused by numerous origins especially due 
to the wind. 

Total range of these projectiles are rather great. Consequently, 
perturbing influences from these origins are important to notice when 
long range trajectories are calculated. Some of the primary origins of 
error are mass unbalance due to impreciseness of construction and from 
shock of launching, misalignment of jet, misalignment of intake air flow, 

a (D 

and so on. 

In order to minimize these effects, it is useful to set the value of 
thrust of the ramjet equal and opposite to that of aerodynamic drag. 

The trajectory is equivalent to projectile flight in a vacuum and is 
known as thrust- equal— drag trajectory. The effect of wind, when the 
direction of wind is the same as that of the trajectory, can be readily 
calculated. In case of a side wind, the projectile will be carried 
laterally. 
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In order to cancel the effect of a side wind which will carry the 
projectile laterally and to restore the planar motion of the trajectory, 
active control can be used. For example, canards can be employed. Due 
to active aerodynamic control, total range will be affected. Active 
control modulates aerodynamic drag. The boundary and initial conditions 
will influence whether total range will be expanded or shortened by the 
effect of canards. 

In this paper, effects due to winds are considered; however, pertur- 
bations of other origins are concentrated in one term in the equations of 
motion as constants. In order to calculate these data by electronic 
computer, flow diagrams for FORTRAN are indicated. 

2. Coordinate Systems 

In order to construct solid body equations, coordinate systems are 

( 2 ) 

considered as follows: 

(1) Ground axis coordinate system (0-xyz). 

The origin is on the center of gravity of the projectile at the 
time of departure of projectile. The x axis is toward the 
direction of horizontal line of shot. The y axis is vertically 
upright, and the z axis is horizontally rectangular in a right 
hand sense. 

(2) Velocity axis coordinate system (0-x^y^z^) . 

The origin is at the center of gravity, CG, of the projectile. 

The x^ axis is along the velocity vector of the projectile, i.e., 
x^ is tangent to the trajectory. The transformation to create 
the (x v ,y v ,z^) coordinate system from the (x,y,z) coordinate 
system is, first, a rotation about the z axis by an angle 0. 

A new coordinate system is formed and is identified as the (x^,y^,.z) 
coordinate system. 
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about 



The velocity vector, v and x^ are located in the x^-z plane. 
Now rotate the (x^,y^,z) coordinate system about the y^ axis 
until the x^ axis is parallel with v. The angle of rotation 
the y^ axis is a. The resulting coordinate system is the 
(x v , y^, z^) system. Refer to Figure 1 for an illustration of 
the axes. 



. y 




Figure 1. Illustration of geometry for the velocity axis coordinate 
system. 
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(3) 



Nonrotating body axis coordinate system (0 - x^y^Zg). 

The origin is at the center of gravity of projectile: the 

Xg axis coincides with the longitudinal axis of synmetry. The 
y „ and z_ axis are transformed as the same manner as that of 

D D 

velocity axis coordinate system by <J> instead of 9, and by £ in- 
stead of a. 

(4) Wind axis coordinate system (0 - x y z )* The system moves with 

J varva va 1 

relative wind speed similar to velocity axis coordinate system. 

The relevant angles are (0 ,0 ) instead of (6,0) with the x 

a* a * va 

axis along the relative velocity vector. 

3* Equations 

Solid body equations with six degrees of freedom are developed by 
the method of coordinate transformation in order to calculate the trajectory 
of the projectile with ramjet propulsion as follows: 



V cos 0 cos 0 


(1) 


v sin 9 cos a 


(2) 


-v sin a 


(3) 


— F 


(4) 


m xv 


1 F 


(5) 


mv cos a yv 




i 
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(6) 


mv zv 


r / cos £ 


(7) 


q 


(8) 


r tan £ 


(9) 


M n /I 
xB x 


(10) 


(1 - I/Dpr + M yB /I 


(11) 


-a - yi) M + m zB /i 


(12) 
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(13) 



xv 



= -mg sin 9 cos O + T cos (<P - 9) - %pv Z SC 



F = -mg cos 9 - Ta + kpv SC a + C cos it 

y v Lt 0 

a 



zv 



= -mg sin 9 sin a - T $+ %)v Z SC T 3 + C sinir 



M = 0 
xB 

2 

M _ = %pv SdC 3 - Tr 8 + M cos it 
ys m t £ 

2 a 

M - = hpv SdC a - Tr a + M sin tt 
zd m t £ 

a 



a 



a = tan 



-i 



cos ?sin(<jP-6) 



coscr cos? cos (<p-Q) + sino sin? 



3 = tan - ^ sine cos? cos(C9-9) - coscr sin^ 
cosP cos? cos (<p-Q) + sine sin? 



(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 



For a projectile with ramjet propulsion, eq. (16) is true usually. At 
launch, p = 0, and there is no initial spin so that p = 0. For case of 
p = 0, equations (11) and (12) become 



q ■ M yB /I (11) ' 

r-M zB /I (12)' 

Furthermore, if thrust has a value near to drag, then dispersion becomes 
a minimum so that the last terms of equations (14), (15), (17) and 
(18) are deleted from these equations. 

Effect of Wind 

When a wind blows, a wind axis coordinate system is considered. 

Using these coordinates, the x axis is along the relative velocity between 
projectile and air. The motion of the projectile can be calculated the 
same way as in the case of no wind. For this case, we use a subscript (a) 
which is applied to the notations of the case of no wind. 

In order to obtain velocities which were observed from the ground, 
we must add distances transported by wind during the time interval to 
values calculated on the wind axis coordinate. The relations between 
velocities on wind axis coordinate and velocity axis coordinate are as 
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a - v cos 9 
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sin 
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a = v sin 9 
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(22) 
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sin 
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= V 


sin a + w 






(23) 


V a 




'a 




z 








From these equations. 


we 


can 


calculate v 

a 


’ a’ a 


knowing 


data without 



wind, v, 9, a and the component of wind relative to the ground w x W y W z * 

In the case of computer calculation by a step-by-step method we 
can use (21), (22), (23) to get initial value of v, 9, O in each suc- 
cessive step. These values of v, 0, a are relative to the ground axis 
coordinates, so that the trajectory is obtained from them. 

4. Effect of Canard 

From mentioned above, it is clear that projectile is not within the 

plane of launching but flies in a different direction from the first 

schedule when wind is blowing. In order to cancel the miss distance, 

we can use a canard giving some suitable trim angle. Trim angle will be 

given on body axis coordinate and whose effect will act as sinusoidal 

form such as V sin y where y is factor angle measured from y axis right 

wise around x^ axis. In order to cancel the effect of sidewise dispersion 

. « 

z = 0 so a = 0 , a. ='0 must be satisfied, 
therefore 



In order to demand 
terms of eq. (15) 



F =0 
zv 



(24) 



(24) we must add a term for the canard on the right 



2 

-mg sin 0 sin CT - TS + %)v SC 0 + C sin TT = v sin y (25) 

ij £ 

a 



There are many factors to give effect to the left side terms of this 
equation, so it is difficult to discuss analytically However, the 
value of the terms are calculated on each step so that we can decide a 
suitable value of canard angle for each step of the calculation. 
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5, Program Flowcharts 

In order to calculate these equations, a computer program will 
be built. (3 lt is illustrated by program flowcharts* This computer program 
consists of a main program and subroutines. 

Inputs which are used for calculation are ready and stored in incre- 
ment address by orders in the program to use for calculation. Data 
calculated are provided printed on sheets and so on. Program flowcharts 
of main program and subroutines are shown on figures, and they are explained 
as follows; 

a. Main Program 

Main program has a job to store and supply data from incre- 
ment address of computer and to make initial condition in shape 
suitable for calculation. It gives these data of subroutines to com- 
pute and to get results of calculation. Data returned from sub- 
routine are printed periodically. Calculation is carried by one 
step of variables. 

These data are accumulated step-by-step until they come 
to some objective point. When the objective point is satisfied, 
the calculation is stopped. The point where the projectile 
touched the ground, point of print and points to change initial 
conditions are detected and ordered to computer. It contains FORMAT 
to print these results, 
b. Subroutine COMPUTE 

In order to calculate eq. (1) ^ (12) there need data of eq. (13) 

^ (20) which are in the state suitable for input. In this subroutine 
these data are prepared using three subroutines of ATOM, FORCE and 
MO MEN. These data are sent to subroutine RKG. There eq. (1) ^ (12) 
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are calculated and results are sent back to this routine. Additional 



data are calculated here. These data are sent back to main program, 
t. Subroutine RKG 

This subroutine has a job to calculate each of the dif- 
ferential equations in a Runga-Kutta-Gill method. Suitable data 
for input are sent from subroutine COMPUTE. This subroutine calls 
differential equations one by one from subroutine DANDO. 

d. Subroutine DANDO 

This subroutine stores each equations whwch composes the tra- 
jectory equations in state of FORTRAN form. 

e . ATOM 

This subroutine provides air temperature and air density for 
the standard state of heights from sea level. Air temperature is 
used to calculate Mach number. 

f. FORCE 

Data of forces are stored in increment address. On this sub- 
routine forces that are arranged suitably to the state of input of 
the instance are selected and prepared for supply. This subroutine 
has a function to make linear interpolation. 

g. MO MEN 

Data of moments are stored in increment address. Moments that 
are arranged suitably to the state of input of the instance are 
selected and prepared to supply. This subroutine has a function 
to make linear interpolation. 
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Program Flow Charts 



(a) Main Program 
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(b) Subroutine COMPUTE 
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(C) Subroutine FKG 
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v ✓ 
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(d) Subroutine DAN DO 
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(e) Subroutine ATOM 
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(f) Subroutine FORCE 
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(g) 



Subroutine MOMEN 
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6. Conclusions 



Solid body equations with 6-degree-of-f reedom are formulated in order 
to calculate the trajectory of a projectile which is flying with the thrust 
of a ramjet. 

In order to calculate by electronic digital computer; program flow 
charts are developed* FORTRAN programs will be made for each subroutine and 
program. Trajectory errors due to the effect of the wind are considered 
mainly. Dispersion due to the wind will be calculated by these programs. 

If there need to cancel the effect of a side wind, trim angle of canard will 
be given step by step in the procedure of calculation. 
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